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Abstract. We compute the 0(1/ Nf) correction to the predominantly gluonic flavour singlet 
twist- 2 anomalous dimension used in polarized deep inelastic scattering. It is consistent with 
known two loop perturbation theory and we determine the three loop contribution at 0(1/Nr). 
The treatment of the e-tensor in the large Nf d-dimensional critical point formalism is also 
discussed. 
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There has been intense experimental activity in recent years into understanding the origin 
of the spin of the nucleons. This is primarily because of the EMC effect, Q, which appeared 
to suggest that the proton spin was due to a sea quark effect with no significant contribution 
from the valence quarks. Consequently there has been theoretical interest in the problem. For 
example, one avenue of research has been to improve theoretical predictions from perturbative 
QCD by going to higher loop orders, 0. A problem related to this concerns the construction 
of the perturbative sector of the operator product expansion used in polarized deep inelastic 
scattering which is important, for example, for current and future HERA experiments. Unlike 
the unpolarized sector, polarized calculations of twist-2 operator dimensions have only recently 
become available at two loops, || ||], as a function of the operator moment, n, or equivalently 
as a function of x, the conjugate variable of the Mellin transform in the DGLAP formalism 
of Dokshitzer, Gribov, Lipatov, Altarelli and Parisi, [||, ||. Although the two loop results of 
H Q represent progress for polarized scattering, building on the one loop results of Q, there is 
clearly a need to push calculations as far as is possible. Indeed for unpolarized scattering the 
relevant twist-2 operator dimensions are available exactly at three loops for the even moments, 
n < 8, (and additionally n = 10 for the flavour non-singlet case, ||). Although the explicit 
n-dependence of the three loop result is the main goal, since it is needed to carry out the full two 
loop renormalization group evolution, some progress has been made in that direction already. 
Recently, the finite parts of the two loop graphs have been determined in |l(| as a function of x. 
These represent the first stage of a full three loop calculation since they will contribute to the 
final value when multiplied by one loop counterterms. Further, some information is available 
from the large Nf expansion where the 1/Nf leading order non-singlet anomalous dimensions 
have been deduced to all orders in the coupling, explicitly as a function of n, [jllj] , with Nf 
the number of quark flavours. (The method used is based on a critical point argument where 
one considers the theory in the neighbourhood of a d-dimensional non-trivial zero of the (5- 
function and computes critical exponents. The initial development was for the O(N) a model 
in p2|. Through the critical renormalization group equation these exponents are related to the 



perturbative renormalization group functions.) The restriction of 11] to compare with the exact 
three loop non-sing let results of g gave agreement where there is overlap. In anticipation of a 
full polarized calculation at three loops the large Nf unpolarized calculations have been extended 
in |l3[| to the polarized case. However, in the calculation of [|l3| only the dimension of the twist-2 
flavour singlet operator which was predominantly fermionic in nature, was considered. The aim 
of this article is to complete the gap in the twist-2 singlet polarized sector by computing the first 
non-trivial large Nf contribution to the flavour singlet operator which is gluonic in character. 
In providing this result we will not only be able to check the recent two loop results but also 
be able to provide the first insight into the n-dependence of the three loop dimension, albeit at 
0(1/Nf). As the calculation of the unpolarized singlet gluonic operators was treated at length in 
we will be brief concerning the technical details of the integration and focus on the salient 
differences which arise in the polarized case. 

We recall that the basic twist-2 flavour singlet polarized operators are, ]?]], 

Of = i n - 1 5^7 5 7 w L" t2 ...D^ip - trace terms 

= li^Se^^^G^D^ ...D^G a ^ a - trace terms (1) 

where if) is the quark field, G£„ = d^A a u - d v A^ + f abc A b ^A c v , A% is the gluon field, f abc are the 
colour group structure functions, S denotes symmetrization over the Lorentz indices and a^uap 
is a totally antisymmetric four dimensional pseudotensor. As has been observed in j/J, since O q 
and O g have the same dimension, they will mix under renormalization in perturbation theory. 
This leads to a mixing matrix of renormalization constants and hence anomalous dimensions 
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which we denote by 

( 7 99 (o) Igqifi) \ 
V 7g 9 (a) 7 99 ( a ) / 
where we use the same definition for the coupling constant, a, as 
Nf dependence of each of the entries will be 



(2) 

To three loops the explicit 



7<?<?(«) 
7<?<z( a ) 

lgg( a ) 



aia + (a 21 N f + a 22 )a 2 + {a u Nf + a 32 N f + d 33 )a 3 + 0{a 4 ) 
ha + (b 21 N f + b 22 )a 2 + (631^/ + 632^/ + b 33 )a 3 + 0(a 4 ) 
dNfd + c 2 N f a 2 + {c 3X Nf 
(dnNf + d\ 2 )a 



f 

c 32 Nf 
,2 



C33ja 



0(a 4 



(d 2 iJV/ + d 22 )a 2 + (d 3 iJV/ + cfeiV) + d 33 )a 3 



0(a 4 



(3) 



where Nf = T(R)Nf and the coefficients aij , 6^-, Cjj and djj will depend on n and the colour group 
Casimirs. This form will be important when we apply the critical point formalism to deduce 
the 0{\/Nf) anomalous dimensions. Moreover, for concreteness and to clarify the conventions 
we use, we recall the one loop polarized results of Ahmed and Ross are, 0, 
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ci 



C 2 (R) 
4(n 



45i(n) - 
- 1)T(R) 



du 



n(n + 1 
4 



2(n + 2)C 2 (R) 
n(n + 1) 



T(R) 



(4) 



n(n + 1) 3 

where the finite sums Si(n) are defined by Si(n) = J2r=i 1 / rl - Although we have indicated the 
form of the mixing matrix in (Q), one could in principle consider a basis where it is diagonal. 
Indeed it turns out that in the 1/iVf formalism of |l3| that that method produces the anomalous 
dimensions of the twist- 2 eigenoperators of (g). This is primarily due to the fact that the 
canonical dimensions of the operators of QCD in perturbation theory are the same, whereas 
computing at the fixed point of QCD in (/-dimensional spacetime, where the large Nf expansion 
is constructed, the canonical dimensions are different. Indeed it is elementary to deduce this 
difference is 0(e) where d = 4 — 2e. Further, it was shown in [13] that the critical exponent of 
the operator Of corresponded to the eigen-anomalous dimension of ([!]) with a predominantly 
fermionic content evaluated at the non-trivial zero of the j3- function. Likewise in [14], it was 
demonstrated that the remaining operator corresponded to the predominantly gluonic operator. 
More concretely, diagonalizing d2|) produces the eigen-anomalous dimensions 
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1 

The critical value of the QCD coupling has been given before in |14[] and is 
3e 
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+ 
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qq 9Q7 
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where we have included the contribution from the recently computed four loop /^-function, [ jD| . 
Therefore substituting (||) in (||) and expanding in powers of 1/Nf, we find that the gluonic 
dimension to third order in e is given by 
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du 4 
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To reproduce this at 0(e 2 ) together with the subsequent terms at 0(1/Nf) one applies the 
critical point formalism to the operator O g which exploits the fixed point equivalence of QCD 
with the non-abelian Thirring model, Jig] , at leading order in 1/Nf. As the method has been 
well documented, jnj, we will comment only on those technical details of the computation which 
are peculiar to O g . Indeed as the calculation virtually parallels that of the unpolarized case we 
refer the interested reader to Q]. There details of the integration procedure for the Feynman 
diagrams have been recorded, including the construction and solution of the recurrence relations 
for deducing the contributions from subintegrals in the graphs. These have been implemented in 



the symbolic manipulation language form, [17], which was used to perform the tedious algebra 



in the integration. The four contributing Feynman graphs for the 0(1/Nf) part of (0) are 
illustrated in fig. 1, where the operator insertion is indicated on an internal gluon line. We have 
used the same Feynman rule for the insertion and conventions as M. 



One technical issue deserves detailed comment. As has been noted elsewhere, [lq] , one 
has to deal with the problem of treating an object, such as c f^vcr p<, which is defined only in 
four dimensions, in a formalism which is manifestly d-dimensional. Such a situation arises 
not only in large Nf calculations but also in explicit perturbative multiloop computations. Its 
resolution in the latter case is to treat products of two e-tensors, where they arise, as the 
determinant of a 4 x 4 matrix of r\p V tensors, p9| . As the latter can have a c?-dimensional 
interpretation one can perform a calculation without reference to the original e-tensors. To 
compensate for the obvious discrepancy and to restore the Ward identity, one introduces an 



additional (finite) renormalization constant, [2C[. This approach was developed for the treatment 



of the renormalization of the singlet axial vector current which has an anomaly, p9, 21] 



similar approach is adopted for 1/A^- calculations, [22], which has already been applied to the 
calculation of the exponent corresponding to the anomalous dimension of the axial vector current, 
|l3|[ . For each graph of the present calculation, we project out the contribution to the exponent 
by multiplying by e Atl/0 -pA f7 p p , where is a null 4- vector used to define the Feynman rules on 
the light cone and p 11 is the external momentum flowing through each diagram. Therefore each 
integral contains the product of two e-tensors and like p0[| we replace these by a determinant of 
tensors and compute in d-dimensions. However, we also need to introduce a compensating 
factor since not only is ep U<J p four dimensional but the vectors A^ and p^ are strictly physical 
objects and cannot be extended to <i-dimensions, [l^]. In the split 7-algebra notation of [113], A^' 
= A' 2 and p^ = p^ where the bar denotes the four dimensional subspace of the d-dimensional 
extended spacetime, which corresponds to the physical spacetime in the limit as d — > 4. To 
summarize we have replaced the correct e-tensor combination given on the right side of 

e w ^AVAV = /(^ Fi pe aW -AVAV (8) 

by the object on the left side which is formally d-dimensional and involves the //-determinant. 
The discrepancy in this procedure is represented by a compensating function of the spacetime 
dimension. We have computed this proportionality function, /(/i), by considering contractions 
of the (i-dimensional object and comparing it with similar contractions of the object restricted 
to 4-dimensional spacetime. We found 

/Ox) = (2/i-3)(/i-l) (9) 

To check that this is a correct procedure we computed the first graph of fig. 1 by using the 
determinant substitution approach and then compared the result with that obtained by direct 
calculation of the integral. Such a computation is possible in this case since the graph is in 
effect a series of chain integrals in the language of |l^]. It turns out that by inspecting each of 
the other graphs a combination which is effectively the same as (||) arises and therefore (|9j) is 
used for each of the graphs. Moreover, this argument was also used to justify the inclusion of a 



4 



similar factor in the large Nf calculation of the anomalous dimension of the singlet axial vector 
current, 13]. The final expression for that case was consistent with three loop perturbation 
theory, 0,0. 

In light of these remarks, we can now quote the final value we obtained for A+ J i(a c ) where 
A + (a c ) = J2iZo ^+,i{ a c)/Nj. After evaluating the contribution from each graph we find 

(n + 2)(n - l)T(n + 2 - (j,)T(n + l)C 2 (R)v? 



A+,i(a c 



+ 



(/x - 2) 2 (n + 1)T(2 - ti)T{pi + n)nT(R) 

2/i(/x-l)Si(n)C 2 (G>ff 
(2/x-l)(^-2)T(i?) 



- [4n 3 n 2 + 4/i 3 n - 8/x 3 



8/iV 



8^ 2 n + 16/i 2 + 



+ 5/m — 9/u. 



n 



n + 2]r(n + 2 - p)r(/i)/iC 2 (G)»7? 



/[8(2/x — 1)0*— 1) 3 (^ + l)r(3 - /i)T(/i + n)nr(i?)] 



[32/i 4 - 4/i 3 n 2 - 4/i 3 n - 120/^ 3 + 16/iV + 16/n 2 n 
+ 160/j 2 - 20/in 2 - 20/m - 89// + 8n 2 + 8n + 18] nC 2 {G)r^ 
/[8(2/x - l)(/i - l) 3 (/i - 2)(n + l)nT(R)] (10) 



where we have used reduce, P3| , to handle tedious amounts of algebra. The quantity 77° is 
defined by 

o _ (2/x-l)(/i-2)r(2/i) 

771 4r2( /U )r( M + i)r(2- M ) 1 j 

There are various checks on the validity of (|io|). First, we find that expanding ([h]) in powers 
of e it agrees with two loop perturbation theory, jj], ||, From one point of view this justifies 
the use of (|8|), since (^) is e-dependent and therefore will affect the form of the 0(e 2 ) terms 
of (|l0|) when it is expanded. Second, if one substitutes n = 1 in (|To|), we find that it vanishes 
identically. This is consistent with the observations of [ 24 , 25 1 , where it was shown that the axial 



anomaly forces the n = 1 gluonic operator to have vanishing anomalous dimension. Therefore 
we are confident that (pi) is correct and can make an observation on the form of the 3-loop 



contribution to the gluonic eigen-anomalous dimension. By expanding (10) to 0(e 3 ) and using 
(0), we find that 

&31CL 64(7n 2 + 7n + 3)(n + 2)(n-l)S 1 (n)C 2 (R) 
«3i + — ; — — 



+ 



in 9(n + l) 3 n 3 

Q4(n + 2)(n-l)S 2 (n)C 2 (R) 
3(n + l) 2 n 2 

4[33n 8 + 132n 7 + 142n 6 - 36n 5 - 263n 4 - 312n 3 
+ 280n 2 + 408n + 144]C 2 (i?)/[27(n + l) 4 n 4 ] 
8(8n 4 + 16n 3 - 19n 2 - 27n + 48)5i(n)C 2 (G) 
27(ra + l) 2 n 2 

2(87n 6 + 261n 5 + 249n 4 + 63n 3 - 76n 2 - 64n - 96)C 2 (G) 
27(n + l) 3 n 3 



(12) 



Evaluating this at n = 1, one obtains a value which, from ([?]), is correctly proportional to the 
three loop coefficient of a c . As an aid to future explicit perturbative calculations similar to 
||], we have evaluated (|l^) for odd n, n < 23, and listed the results in table 1. Unlike in the 
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unpolarized calculations, here it is the odd moments which are important since the original 
operator does not exist when n is even as can be seen, for example, from the Feynman rules, |J. 
A further remark on ( |l2| ) concerns the status of its C2(G) sector. In the unpolarized case, [|l4j], 
it was noted that at 0(1/Nf) 631 will only involve contributions proportional to C2(R) and not 
02(G). Assuming the situation is similar here, this implies that the terms involving 02(G) in 
(|l~2| ) are in fact the correct n-dependent part of the coefficient set ds± of the original perturbative 
operator. Finally, several general remarks on the overall structure are in order. First, the full 
expression is much simpler in form than that of the unpolarized case, Indeed this was also 
observed in the fermionic singlet sector as well. Second, in the evaluation of each integral an 
extra n-dependent structure was present which was proportional to r(u — l)T(n)/T(fi + n). This 



is similar to the unpolarized situation, [14], and like in that calculation, it also cancels in the 
final summation. 

To conclude with we note that all the 1/Nf leading order terms are now available for the 
anomalous dimensions of the twist-2 operators used in polarized deep inelastic scattering. The 
next stage in this 1/iVy programme would be to consider the unpolarized flavour non-singlet 
sector and compute the 0(1/Nj?) correction to the dimension of the twist-2 operator as a function 
of n. 
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Table 1. Coefficients of \dz\ + b^ici/dn] as a function of moment. 
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Fig. 1. Leading order diagrams for A + (a c ). 
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